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We  projxDse  to  discuss  the  problem  of  the  stability  of  a  cylindrical 
panel  resting  on  a  rectangular  boundary,  under  the  effect  of  a  normal  im¬ 
pulse  which  is  characterized  by  a  rapid  increase  in  load  to  a  certain  magni¬ 
tude  and  subsequent  decrease  according  to  the  exponential  law.  Our  main  ob¬ 
jective  is  to  construct  a  range  of  parameters  which  characterizes  this  im¬ 
pulse  and  in  the  presence  of  which  snap-through  buckling  of  the  shell  does 
not  occur  (l.e.,  stability  "in  the  large").  The  effect  of  the  initial  com¬ 
pression  and  damping  decrements  on  the  boundary  of  the  stability  region,  as 
well  aa  on  the  maximum  deflections  attainable  under  the  effect  of  the  im¬ 
pulse  are  investigated.^ 

Partial  differential  equations  of  the  nonlinear  shell  theory  are,  by 
applying  the  Papkovich-Oalerkin  method,  reduced  to  nonlinear  ordinary  dif¬ 
ferential  equations  and  solved  on  continuously-acting  analog  computers  (in 
this  work  the  MN-7  computer  was  used). 

1.  Basic  Equations.  We  shall  use  the  basic  hypotheses  of  the  non¬ 
linear  theory  of  thin  elastic  shells  G;  C  assuming  that  elongations  and 
shears  are  small  in  comparlscxi  with  unity  wd  that  the  deflections  of  the 
shell  are  compEurable  with  its  thickness,  though  small  in  comparison  with  its 
other  dimensions.  We  shall  also  assume  that  the  Kirchhoff-Love  hypothesis 
remains  valid  in  the  case  of  finite  deflections  of  l.he  shell  as  well. 

Under  the  assvnptlons  made,  th(^  deflection  equations  of  the  shell 
with  principal  curvatures  k]^  and  k2  have  the  following  form: 


One  of  the  authors  presented  this  paper  at  the  conference  (3)  on  elastic 
oscillations  at  the  Institute  of  Mechanical  Engineering  of  the  Acadeiqy  of 
Sciences  of  the  Latvian  S.8.R.  (Riga,  June  19^)* 
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Here  Wo(x,y)  Is  the  initial  deflection^  v(x,y,t)  is  the  complementary  de¬ 
flection,  X  and  y  are  the  coordinates  computed  along  the  lines  of  curva¬ 
ture  in  the  middle  surface,  D  is  the  flexural  rigidity,  E  is  the  modulus 
of  elasticity,  h  is  the  thickness  of  the  shell,  *(x,y,t)  is  a  stress  func¬ 
tion  of  the  middle  surface  stress -resultants  related  to  a  unit  length. 
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The  symbol  L(A,B)  denotes  a  bilinear  operator 


L(A,B) 


a^A  a®B  a^A  A  „  a^A  a^B 
a7a7  ^ax^ '  Sizars? 


If  the  shell  oscillates,  then  the  normal  load  component  q  vill  be 
computed  using  the  formula 


q(x#y/t)  -  -  ph 


at" 


2phe  ^  4-  qp(x,y,t) 


(1.3) 


where  p  is  the  density  of  material  of  -che  shell,  €  is  the  damping  coeffi¬ 
cient,  and  qo  is  the  external  load. 

We  note  that  along  with  normal  inertia  forces,  the  tangential  in¬ 
ertia  forces  should  be  considered  in  the  equation  as  well.  However,  in  the 
case  where  loading  is  not  too  fast,  ,  if  the  time  for  the  e.xt.emal  forces 
to  vary  by  an.  appreciable  magnitude  is  suff  ij:ient.l.y  long  in  c.omparison  to  the 
longest  period  of  tangential  natural  oscillat-ions  though  commensurable  with 
the  longest  period  of  lateral  oscillations,  then  tangential  forces  of  inertia 
can  be  disregarded..  An  analogous  assumption  is  usually  introduced  in  the  in¬ 
vestigation  of  natviral  and  forced  oscillations  of  plates  and  shallow 
shells 
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To  solve  the  system  of  equations  (l.l)  for  a  case  of  normal  load 
■«v  defined  by  formula  (I.3)  we  shall  use  the  Papkovlch-Galerkln  method. 

We  shall  represent  the  initial  and  complementary  deflections  in 
the  form  of  series  of  functions  ♦jCx^y)  satisfying  the  boundary  condltlone 


=  Y.  Cj^Vj(x,y) 

J— 1 

flO 

w(x,y,t)  =  Y  ^.(^^.(x^y)  (1.4) 

J=1  •’ 


•<  * 


Here,  5(t)  are  the  unknown  time  functions. 

Functions  ♦j(x,y)  are  selected  in  such  a  way  as  to  satisfactorily 
describe  the  expected  mode  shape  of  the  oscillations  of  the  shell;  in  the 
following  we  shall  assume  that  the  mode  shape  of  the  small  natural  oscilla¬ 
tions  of  the  shell  are  such  functions.  At  first,  we  shall  substitute  the 
series  (1.4)  into  the  second  equation  (l.l)  and  solve  it;  we  shall  then 
satisfy  boimdary  conditions  of  the  stress  function  ®(x,y,t)  and  substitute 
it  into  the  first  equation.  Finally,  we  will  apply  the  Galerkln  variaticmal 
method,  and  the  problem  will  be  reduced  to  a  system  of  ordinary  differential 
equations 


.t)  = 
m 


(o) 
1  ' 


J  “  1>2 


, . . . 


(1.5) 


Here  €j  denotes  damping  coefficients  which,  in  contrast  to  formu¬ 
la  (l.3)>  are  ass\jmed  to  be  distinct,  and  Wj  denotes  the  frequencies  of 
small  natural  oscillations.  Nonlinear  terms  accounting  for  external  load  are 
included  in  functions  f  and  q  . 

J  J 

In  the  case  of  a  shell  with  edges  freely  supported  on  a  contour 
with  sides  a  and  b  (Figure  l),  the  solution  can  be  written  in  the  form 


-  T.  t  sin  sin  2^ 

°  m=l  n=l  a  b 


w(x,y,t) 


The  munber  of  the  retained 
terms  of  the  series  detennlnes  the 
order  of  the  differential  system 
and  thus  is  ccjntroiled  by  the  capacity 
of  the  analog  computer.  The  MN-7  com¬ 
puter  can  solve  nonlinear  equations  of 
the  second  order  •  accordingny  we  write 
the  solution  in  the  form 

w  (x,y)  =  5  sin  —  sin 
w(x,y,t)  =  C(t)  sin  ~  sin  (l.6) 


Substituting  expression  (l.6)  into  the  second  of  the  equations 
(l.l)  we  obtain  equations  of  the  following  type: 


=  F(x,y,t) 


(here,  time  t  is  considered  as  a  parameter). 

Let  us  consider  the  case  of  a  cylindrical  shell  (k-j^  =  0,k2  =  1/R) 
with  the  following  boundeury  conditions:  edges  x  =  0  and  x  »  a  are  under 
^  compressive  stresses  =  -p,  while  edges  y  =»^  0  and  y  *=  b  are  restrain¬ 

ed  from  displacements.  The  tangential  stresses  on  the  edges  will  be  con¬ 
sidered  equal  to  zero.  If  one  takes  for  4  the  expression 
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(1.7) 


then  these  boxmdary  conditions  will  be  satisfied  "in  the  mean. " 

Here,  p^  and  Py  are  constants  which  can  be  deteralned  from  the 

conditions 


— 7T  dy 

dx^' 
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dx  dy  =  0 


(1.8) 
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The  meaning  of  the  fi.rst  condition  is  explicit.  The  second  con¬ 
dition  requires  that  the  displacement  of  edges  y  -  0  and  y  =  b  in  the  di¬ 
rection  of  y  axis  is,  "on  the  average, "  equal  to  zero. 

Substituting  expressions  (1.6)  and  (1.7)  into  conditions  (1.8)  we 
obtain  constants  pj^  and  Py,  after  which  the  function  (l.?)  can  be  com¬ 
pletely  determined. 

We  then  substitute  the  expression  for  this  fvinctlon  into  the  first 
of  equations  (l.l).  In  order  to  reduce  this  equation  to  an  ordinary  differ¬ 
ential  equation  we  substitute  into  it  the  expressions  (1.6)  and  we  apply  the 
Galerkin  variational  method.  After  the  usual  computation,  which  is  not  shown 
here,  we  obtain  the  following  equation: 


T.  ^ 


Vu)2  +  16 


<lo*(t) 

2 

nph 


(1.9) 


Here,  w  is  the  frequency  of  small  natural  oscillations  of  an  un¬ 
loaded  shell,  5(t)  is  a  function  characterizing  the  response,  f(C>C  )  la  a 
nonlinear  function  which  is  a  part  of  the  representation  of  initial  md  sub¬ 
sequent  deflections,  pg  is  a  magnitude  having  linear  stress  dimension,  and 
qo  is  a  certain  constant  characterizing  the  magnitude  of  normal  load  (for 
example,  its  value  for  t  -  O). 
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(1.13) 
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variables: 


For  what  follows^  we  shall  Introduce  the  followlog  dimensionless 


i.r, 


2n€  « 

—  ■  -  m  0» 
W  ^ 


l6q^ 


2  2  2 
*  w  ph 


(1.15) 


The  meaning  of  the  first  three  variables  is  obvious.  The  next  to 
the  last  variable  is  a  decrement  of  the  small  natiiral  oscillations  and  the 
last  variable  represents  a  deflection  of  the  center  of  the  shell  relative  to 
its  thickness  under  the  effect  of  a  load  qo  applied  statically.  The  equa¬ 
tion  (1.9)  takes  the  form  , 


aC“ 


+ 


75 


^  C,  -  ^  +  C'q>(T)  (1.16) 

*^o  o 


Here,  and  in  that  which  follows,  the  bar  over  ^  will  be  omitted. 

2.  Results  of  the  Solution  Obtained  by  Means  of  an  Electronic  Analog 
Computer.  The  equation  (I.16)  will  be  solved  for  the  case  where  the  impulse 
has  an  exponential  form,  i.e.,  for  the  case,  where 


-CT 


<P(t)  =  e 


for  T  k  0(wt  =  t) 


(2.1) 
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Here  c  is  a  constant,  characterlzinf^  the  rate  of  decrease  of  load 
Intensity.  For  example,  If  cn  «•  0.5;  this  raeans  that  intensity  of  the  load 
decreased  hy  a  factor  e  within  one  period  of  the  free  linear  oscillatJ.ons. 
According  to  formulae  (2.l),  the  l0£id  q  increases  instantly  from  0  to 
at  the  instant  of  time  t  =  0.  In  addition,  it  is  necessary  to  take  into 
account  both  normal  and  tangential  inertia  forces.  A  contradiction  will  be 
eliminated  if  one  considers  that  the  expression  describes  an  idealized  pic¬ 
ture  in  comparison  with  actual  conditions  when  pressure  Increases  rapidly, 
yet  with  a  finite  velocity  (Figure  2), 


We  shall  not  discuss  here  the 
effect  of  tangential  Inertia  forces  and 
the  effects  of  elastic  wave  prcq)agation 
occurring  with  large  deflections  of  the 
shell. 


If  an  elastic  shell  is  sub¬ 
jected  to  impulsive  loads,  the  follow¬ 
ing  problems  originate: 

a)  to  find  the  conditions 
vinder  which  buckling  of  the  shell  would 
not  occur. 


Figure  2 


b)  if  buckling  is  Imixissible,  then  to  find  the  greatest  deflec¬ 
tions  and  stresses  originating  with  an  Impact, 


c)  if  buckling  occurs  all  the  same,  then  to  find  the  greatest  de¬ 
flections  and  stresses  which  occur  after  buckling. 

Computations  were  made  for  a  cylindrical  pcmel  resting  on  a  square 
contour  ,  (n  =  l),  with  p  =  0.3  and  k  =  12.  In  this  case,  formulae  (l.l4) 
yield 


a  =  0.162,  p  =  0.485Cq  -  0.690,  7  =  0.3245^  -  0.905Cjj 

When  impacted,  the  shell's  behavior  depends  essentially  on  the  mag¬ 
nitude  of  the  initial  load  p.  In  order  to  investigate  this  problem  we  shall 
first  discuss  equation  (I.16)  for  the  case  of  a  static  load.  For  an  ideal 
shell  (Cq  ®  equation  takes  the  form 

ji  -  ^jc  +  ac^  +  =  5'  -  ^  V  (2.2) 

Dependence  of  the  dimensionless  deflection  ^  on  for  various 
ratios  v/vo  are  shown  in  Figure  3*  H;  follows  from  this  graph  that  for 
sufficiently  small  values  of  p  (approximately  for  p  <  0.3Pq)  and  for 
C'  =  0,  there  exists  one  elastic  equilibrivan  state  of  the  shell.  For 
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P  >  0'3Po  =  0  there  iii;e  thi’ee  i;t|ull:I Im  Iiot  atateHj  one  of  them  cor- 

reeponds  to  the  negative  root  of  equation  (2.  I,  e, ,  -  to  a  deflection  of  the 
shell  with  an  increase  of  its  curvature;  the  other  two  corrtispond  to  the 
positive  roots  of  this  equation. 

The  intermediate  root  corresponds  obviously  to  the  unstable  form 
of  equilibrium.  We  note  that  the  upper  critical  pressure  does  not,  in  the 
usual  sense  of  the  word,  exist  in  this  cases,  since  as  a  result  of  Poisson's 
effect  the  shell  deflections  resulting  from  curvature  Increase  are  present 
even  with  small  values  of  p. 

For  5^  =  0,  the  roots  of  equation  (2.2)  become,  with  T  0,  the 
singular  points  of  equation  (l. l6).  Thus,  for  a  square  panel  and  with  p  < 
0.3Po  there  is  one  slnguleu:  point  (stable  focus),  while  with  p  >  O.SPq 
there  are  three  singular  points  (two  stable  foci  and  a  saddle).  With  p  < 
0.3Poj  and  regardless  of  magnitude  of  5^,  there  is  only  one  type  of  motion 
possible,  namely,  oscillations  with  a  gradually  decreasing  amplitude  around 
a  unique  stable  focus,  l.e.,  oscillations  without  buckling.  With  p  >  0.3Pq, 
the  following  three  types  of  motion  ore  possibles  l)  for  small  5'  —  oscil¬ 
lations  without  buckling  around  a  focus  nearest  to  the  origin  of  coordinates; 
2)  for  large  --  buckling  with  subsequent  oscillations  aroimd  the  second 
focus;  3)  for  very  large  --  buckling  with  subsequent  oscillations  along 
phase  trajectories  embracing  all  three  slnsjular  points.  In  the  iMt  case, 
as  t  >  •,  generally  speaking,  convergence  to  any  of  the  two  stable  foci  is 
possible. 


Figure  3 
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Figure  k  shows  behavior  of  the  shell  with  =0,  p  ■  0. 5Po>  6  ■ 
O.O^n,  and  c«  =2.0. 

If  5^  »  0.68,  buckling  does  not  occur;  with  =  0.72,  there  is 
buckling  with  oscillations  aromid  the  second  focus;  with  “  1*50,  oscil¬ 
lations  occur  with  a  change  of  deflection  sign. 

The  phase  plane  diagram  of  this  case  is  shown  in  Figure  5  (after 
several  oscillation  cycles,  damping  with  a  higher  than  critical  magnitude  was 
added;  this  was  done  in  order  to  accelerate  the  approach  of  the  system  to  a 
stable  singular  point). 

It  is  of  interest  to  construct 
boundaries  of  a  c,^'  domain  of  param¬ 
eters  with  such  a  characteristic  that 
when  the  shell  is  under  the  effect  of  an 
impulse  having  parameters  lying  in  this 
domain,  buckling  of  the  shell  will  not 
occur.  Figure  6  shows  the  boundaries  of 
such  stability  regions  for  an  ideal 
shell  for  various  magnitudes  of  a  longi¬ 
tudinal  compressive  stress.  It  is  clear 
that  stability  regions  are  located  on 
the  side  of  smaller  values  of  5^. 


Figure  5 
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Figure  7 


Figure  8 


It  would  be  of  Interest  to  learn  the  dependence  of  the  critical  im¬ 
pulse  on  other  parameters  of  the  problem.  It  is  easy  to  see  that  a  dimen¬ 
sionless  impulse  I  can  be  determined  from  the  condition 


I 


The  graph  shown  in  Figure  7  vas  plotted  according  to  this  formula. 

Figure  8  shows  the  effect  of  damping  decrement  6  on  the  boundar¬ 
ies  of  the  stability  regions.'  As  was  expected,  this  dependence  turns  out  to 
be  weak,  particularly  when  damping  is  small.  Thus,  when  damping  increases 
from  6.  a  0.01  to  8  =  0. 50,  the  critical  value  of  parameter  increases 
no  more  than 

Figure  9  and  Figure  10  show  the  dependence  of  the  greatest  deflec¬ 
tion  Cmax  of  "the  shell  on  the  intensity  of  a  normal  load  ^  for  values 
p  =  0.2pq  and  p  =  0.5Po'  ^n  the  second  case,  during  transition  of 
through  the  critical  value,  the  maximum  deflections  increase  sheurplyj  how¬ 
ever,  with  sufficiently  strong  impulses  they  differ,  in  both  cases,  only 
slightly  from  each  other.  On  the  other  hand,  the  difference  in  behavior  of 
the  shell  remains  significant  with  t  >  ».  In  this  case,  for  p  <  0. 3Pq, 
and  after  damping  of  oscillations,  the  shell  returns  to  the  initial  state  of 
equilibrium;  however,  for  p  =  0.3Po>  'the  most  probable  form  of  equilibrium 
of  the  shell  after  damping  of  oscillations  appears  to  be  the  buckled  state. 
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Figure  10 


From  the  above  a  conclusion  can  be  drawn  that  an  elastic  shell  can 
withstand,  without  buckling.  Impulsive  loads  exceeding  several  times  the  cri¬ 
tical  static  load.  Let  us,  for  example,  consider  a  case  where  p  =  0.3Po* 

A  static  load  which  originates  buckling  amoimts  (in  dimensionless  magnitudes) 

to  =  0.33-  If  the  load  Is  Impulsive  and  If  cit  =  1.0  (during  the  i>erlod 

of  natural  small  oscillations  of  the  shell  the  load  Intensity  decreases  by  a 
factor  of  e)  then,  as  It  Is  shown  In  Figure  6,  the  dynamic  load  causing 

buckling  Is  =  0.73*  Thus,  the  critical  dynamic  load  Is  2.2  times  greater 

than  the  dynamic  one.  For  loads  whose  Intensity  further  decreases  with  time, 
this  difference  will  be  greater.  For  example,  for  on  ■  ^.0,  the  dynamic 
buckling  load  amounts  to  «  1.28  which  exceeds  by  3.3  times  the  corre¬ 
sponding  static  load. 
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